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ABSTRACT 


An undulatoi is a periodic magnetic field device that is an integral part of a 
Free Electron Laser (FEL). The key to FEL performance is the undulator design. 
This thesis models the undulator magnetic fields by using the Biot-Savart law, a 
current element integration technique. Applying this technique to different 
winding schemes for a bifilar helical compact undulator, the field structure is 
computed for positions outside the entrance to the undulator. These stray fields 
can have a focusing or defocusing effect on the incoming electron beam. This 
disturbs the critical matching of the electron velocity vectors with the co¬ 
propagating laser radiation. The matching of velocity vectors controls the bunching 
of the electrons. This bunching action controls the gain of the FEL, and thus, the 
ultimate performance of the entire device. By choosing a better design for the 
undulator, this unwanted effect is reduced. 

The field structure found from integration determines the electron trajectories. 
The best design for the undulator is then determined by the imposed input 
conditions on the electron beam for entry into the cavity. An additional benefit of 
this technique is a potential application to finding coil winding tolerances for 
undulator construction. 
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I. INTRODUCTION 


One of the designs proposed for a Directed Energy Weapon capable of 
delivering large amounts of energy at a considerable distance at the speed of light 
is the Free Electron Laser (FEL). The FEL uses a relativistic electron beam as a 
power source coupled to the radiation inside of an undulator. The FEL has many 
of the same attributes as ordinary lasers, such as coherence, high rate of fire and 
high energy density [1]. The FEL additionally has the ability to be tuned over a 
large range of wavelengths which increases its flexibility for diverse laboratory 
applications, and as a speed of light weapon, can be tuned for maximum 
atmospheric propagation. 

The key to the performance of the FEL is the undulator design. The 
undulator bends an incoming relativistic electron beam back and forth resulting in 
the emission of photons in the forward direction. If the undulator is placed inside 
a laser resonance cavity, and the dimensions are adjusted appropriately, a FEL is 
made which has the capability to tune the output radiation. The undulator defines 
the magnetic field structure, and thus controls the electron trajectories. Therefore, 
the exact position and velocity of the individual electrons is dependent on the 
design of the undulator. The electron positions and velocities in the co-propagating 
laser radiation are critical to the evolution of electron bunching. This electron 
bunching is responsible for the high energy transfer (gain) between the electron 
beam and the fight. This gain determines important aspects of the performance of 
the FEL. Since the undulator controls the gain achieved, the chosen design of the 
undulator will limit the ultimate performance of the FEL. 

This thesis models the undulator magnetic fields using the Biot-Savart law. 
Using this current element integration technique, we apply it to different winding 
schemes for a bifilar compact undulator. The field structure can then be computed 
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for positions outside the entrance tc the undulator. These stray fields can have a 
focusing or defocusing effect on the incoming electron beam. The matching of the 
electron velocity vectors with the co-propagating laser radiation is essential to the 
energy transfer from the electrons to the light. The energy transfer between them 
is called coupling. This coupling controls the gain of the device. The loss in 
coupling from the focusing or defocusing effect results in lower gain. By choosing a 
better design for the undulator, this unwanted effect is reduced. 

The field structure found from integration is used to determine the electron 
trajectories. The best design for the undulator is determined by the imposed input 
conditions on the electron beam for entry into the cavity. An additional benefit of 
this technique is in showing coil winding tolerances. 

Several new simulation programs have been developed to calculate and 
graphically display the undulator magnetic fields. These calculations were used in 
new simulations for calculating the electron trajectories. This work resulted in two 
talks; one at the Twelfth International Free Electron Laser Conference in Paris, 
France and the other at the International Conference on Lasers ’90 in San Diego, 
California. The talk in Paris has been submitted for publication in Nuclear 
Instruments and Methods in Physical Science, (1991). Additionally, a spread-sheet 
was started to display input and output parameters for comparison of different 
type FELs. The spread-sheet would be useful for a project manager in the Navy to 
determine optimum design parameters as a second check on proposals. This 
project will be completed by a following thesis student. 
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II. BACKGROUND 


A. THEORY OF FREE ELECTRON LASERS 

The Free Electron Laser has many of the same attributes as a regular laser. 
Some of these attributes are "...remarkable directionality, spectral purity, and 
intensity” [2], An ordinary laser is composed of several components (Figure 1). 



Figure 1 . 

Elements of a typical laser from Ref. [2]. 


The essential elements of a laser are: (i) a lasing medium; (ii) a pumping process; 
and (iii) suitable optical feedback elements that allow a beam of radiation to either 
pass once through the laser medium (as in a laser amplifier) or bounce back and 
forth repeatedly through the laser medium (as in a laser oscillator) [2]. 
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The Free Electron Laser is constructed differently (Figure 2) and uses 
different physics to evolve the light amplification. Figure 2 shows the overall 
components and the relationship between the electron trajectory and the co¬ 
propagating light. The FEL lasing medium is a relativistic electron beam. The 
FEL pumping process is from the energy exchange between the electron beam and 
co-propagating light in a vacuum. Mirrors form the optical feedback system shown 
for an oscillator configuration. The amplifier configuration does not have mirrors. 
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INCOMING 

Electrons 


outgoing 

electrons 




Figure 2 . 

Schematic showing basic FEL setup. Also shown is the positional (phase) 
relationship between the light wave and the electron trajectory from Ref. [3]. 
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One of the attributes that the FEL possesses is the capability to "tune" the output 
radiation over a range of wavelengths by changing the electron beam or undulator 
parameters. Another attribute is the absence of extraneous material in the lasing 
medium, allowing the possibility for increasing the power density to weapons grade 
potential. Thus, the FEL is a candidate for a Directed Energy Weapon, and its 
tunability gives it flexibility in optimizing atmospheric propagation [1], 

Given the description of an FEL made earlier, we need to understand the 
physics of the device. Upon entering the undulator, the relativistic electrons are 
"wiggled" back and forth by alternating magnetic fields. From electrodynamic 
theory, we know that a relativistic charged particle that is accelerating produces 
electromagnetic radiation in the forward direction. The radiation is emitted in a 
cone centered around the particle velocity that resembles a "flashlight" (Figure 3). 



Figure 3 . 

Relativistic emission of radiation by an accelerating charged particle. 

Since the width of cone is wide compared to the cavity dimensions, most of the 
radiation is lost and only 10% is retained. Subsequent electrons will radiate into 
the cavity also. Even though there is only about one photon radiated per pass per 
electron, the electron beam density is approximately 10 12 electrons/cm 3 so that 







there is a large photon density. Due to the resulting large photon density, 
spontaneous emission will rapidly be dominated by stimulated emission leading to 
coherence. Spontaneous emission is responsible for the self-starting capability of 
the FEL. As spontaneous emission is replaced by stimulated emission, the 
conditions are set to start matching the electron velocity vectors with the co¬ 
propagating radiation. This matching leads to energy transfer between the 
electrons and the light and is known as coupling. The coupling is responsible for 
the gain of the device. The high energy transfer and efficiency of the FEJ comes 
from the coupling of the free electrons to the coherent radiation set up in the 
cavity. 

To understand this coupling, we need to develop a few equations. Assume the 
radiation is established with a particular wavelength (Figure 4). The velocity of 
the electron is 0* = ]3c where c is the speed of light. Since the electrons are 
relativistic, the transverse velocities are small with respect to the velocity in the z 
direction, |]?|= p z = v z lc . The period of the transverse electron motion 
corresponds to an undulator wavelength X 0 . The condition of resonance occurs 
when the speed difference between the electron and the light is such that one 
wavelength of light, X, passes over the electron as the electron travels through one 
undulator wavelength, X 0 . The electron travels through the undulator period of 
Figure 4 in the time interval A/ = X w /(P 2 c). Considering the difference in speeds, 
the condition of resonance describes the wavelength of light after the electron has 
traveled through an undulator period. At resonance, we find 

X = cAJ-p r cA* = c At-X 0 = (l-p z )VP z • (1) 


For relativistic electrons with small transverse velocity, the Lorentz factor is 



but p 0 = P 2 = — so, p z 
c 


(i - y ~ 2 ) 1/2 


( 2 ) 
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Using the binomial theorem (1 - x) n = 1 - nx ■>. (n/2Xn-l)x 2 - ■ • ■ , and 
neglecting second-order and higher terms, we arrive at 


K- 


A 


fc = 1- 

—H 


2Y 2 


(3) 



S-iart vndulaior 

pe r icd 



period 


e/tc^ron Spe.eA radiation spe&d 

Figure 4 . 

Schematic showing relative positions of the electron and light 
wave through an undulator period ^. 


Using (3), we have 



(4) 


This is the relationship between the undulator wavelength, X 0 , the electron 
energy, y, and the wavelength of light in the cavity, X. If either y, or X 0 , are 
changed, the wavelength of the cavity radiation, X, is changed. This shows the 
tunability of the Free Electron Laser. 

The preceding discussion shows some of the relationships between different 
parameters in the undulator. Figure 5 shows a planar representation of the 
relationship between the electron trajectory and the undulator field. We want to 
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find the equations of motion for the electron so we need to develop expressions for 
the forces and fields affecting it. 


x 



Figure 5 . 

Simplified planar representation of one undulator 
period and resultant electron trajectory. 


The charged electron feels a force given by 

F = q(F + 0* ) , (5) 

where 

v * is the electron velocity, 

F is the electric field felt by the electron, 
q is the charge of the electron, and 
F is the magnetic field felt by the electron. 
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The only fields felt by the electron are due to the cavity radiation, £ r , and the 
total magnetic field which consists of the undulator magnetic field, P m and the 
radiation magnetic field, P r . Therefore, 

P - -j-( y mv* ) = fic ) = q ( P T + f?c x £ ) 

at at 

= + $ X( - B r*+Zli ■ ( 6 ) 

Now, using the definition of momentum gives 

F = f - mc f, m ■ 

but the electron energy is E = ymc 2 . The energy change is cLE !dt = P ■ so 

^ = mc2 iu * qt - e ' = q e r ]tc 

~ ( M,I ■ (7) 

at me 

We now need expressions for the electric and magnetic fields. We want to 
relate the equations to our undulator design, so we choose helical polarization for 
the electric field to match the helical magnetic field undulator design selected. 
From basic optics, the generic forms for planar electric and magnetic radiation 
fields are [4], 


P x (z,t) = 

E x [ cos (kz - <M ) , 0 , 0 ] , 

(8) 

Py(Z,t ) = 

-E y [ 0 , sin (kz - tot ) , 0 ] , 

(9) 


where 

^ , 

c c 

and a) = the angular frequency of the wave. 

Circular polarization can be achieved by combining the planar forms with a spatial 
offset between them called a phase difference, e. The waves move in the positive z 
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direction and combine to form a rotating E vector that describes a helix centered 
on the z axis as shown in Figure 6. If we put this in a combined form, 

£ = £[cos(&z - a it) , -sin (kz - cot) , 0 ] , (10) 

it is easy to see that the magnitude of the electric field vector is constant, but the 
direction of the vector is a function of z and t. 


Y 




Construction of the total electric field vector ( £? ) 
from the component forms. The rotating £* describes 
a circular helix along the z axis from Ref. [4 J. 


If we set \y = (kz - cot), we can simplify the electric field vector to 

E r = E x + E? y = E[ cos(\j/) , —sin (\j/> , 0 ] . (11) 
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Now that we know the form of the electric field, the magnetic field form can be 
deduced. To find 2? r , we use the required conditions for an electromagnetic wave 
that the electric and magnetic field vectors are perpendicular and their cross 
product points in the direction of the direction of propagation. Thus, we can set up 
the following schematic representation and find the orientation and thus the 
solution of the magnetic field vector (Figure 7). 



Figure 7 . 

Schematic showing the orientation of the magnetic field vector, 
, given the electric field vector, F? r . 


The equation for the magnetic field in Figure 7 was arrived at by using the form of 
the electric field and translating it to fit the conditions we just described. From 
geometry, we know that sin(y) = cos(tc/2 - \y) and cos(y) = sin(7b2 - \p). Thus, 
the simplified form for the magnetic field becomes 

S r = — [ sin(vy) , cos(iy) , 0 ] . (12) 

c 
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The helical magnetic field has the form, 


S m = B m [ cos(k Q z ), sin(& 0 z) , 0 ] , (13) 

where k 0 = 2 KlX 0 . 

By substituting equations (11-13) for the fields into our equations of motion, we can 
arrive at the components of the Lorentz Force equation: 


~( yp, ) = —— £(l-p 2 Xcosy , -sin v , 0) + p 2 c£ m (-sin£ 0 z , cos k 0 z , 0) , (14) 
at me L 

-T"( YPz ) = —— E( PjCosy-(Lsiny ) + cP m ( P 2 sin& 0 z - (3 y cosk 0 z ) , (15) 
at me L 


= ——[ P*cosy - Py siny ] 
at me J 


In the equations above, the large scale motion is $ = P z f. The next largest scale 
motion is the transverse "wiggling" motion, P^. For relativistic electrons, 
P 2 —» 1, so ( 1 - p 2 ) -» 0, and we eliminate the first term of (14) above. 

By integrating (14), the transverse velocity ratio is found to be 


= -—[ cos(A 0 z) , sin(fc 0 z) , 0 ] , (17) 

Y 

where K = eB m X 0 /(2nmc). We insert (17) in the relation for the energy change, 
(16), to get 



eKE , Y , 

- cos( £ ) , 

mey 


(18) 


where C, = (k+k 0 )z-(Dt. 

If we take a closer look at the electron phase, C, we see that the important 
factor is the position of the electron. At the beginning of the undulator, time t = 0, 
£(0) = (^ = (k + k 0 )z 0 = 2jez 0 /X, since k ^>k 0 . Thus, the electron phase, £, is a 
"microscopic" variable in that it scales with the optical wavelength. We can relate 
the value of £ to energy exchange between the electron and the light wave. We can 
see that the physics of the energy exchange is proportional to cos( C, ), and 
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-1 < cos( £ ) < 1. The initial electron phases, are random with respect to the 
light wave, and thus, form a uniform distribution of phases. Therefore, half of the 
electron phases, ^, are such that y > 0, and half are such that y < 0. We need a 
better way to express the energy exchange, because it looks like we’re not going to 
get any energy exchange from (18)! Lets try relating y and £ to complete the 
feedback loop. From (2), we have 


\ = 1 _ ( $2 + ft. ^ = 1 ~ P, 2 ~ P l • 

but using (17), we find that p 2 = K 2 ly 2 Therefore, p 2 = 1 - (1 + K 2 )ly 2 . By 
taking the time derivative and some algebraic manipulation, we can find that for 

y> 1, 


y = f lA 

1 + K 2 (k + k Q )C 


f yl 

1 +K 2 kc 


since k ^ k Q 


(19) 


To get this in a better form, it can be shown that by using (3), and the relation 
P 2 = 1 - (1 + K 2 )lf, 


X = -\(1 +K 2 ) , 

2y2 


( 20 ) 


which is a more accurate expression describing FEL resonance than (4)! By 
inserting (20) into (19), we find that 


Y 


eKE 
y me 


cos ( C ) 


2u 0 


Substituting into (18), we have 



2 eKE 0 ) o 
•fmc 


cos( ^ ) 


( 21 ) 


This expression can be made more meaningful by defining a dimensionless time t, 
where x = $ 0 ct/L = ct/L. The dimensionless time has the values 0<x<l 
along the undulator length L. We define t because the total time spent in the 
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undulator is about Lie = 10 -8 seconds for a typical undulator. Using the 
definition of x above, we have 


dx 2 = —-dt 2 and ( ) = 

L z 


d 2 ( ) 
dx 2 


Expression (21) can be simplified to the pendulum equation, 


C = |a |cos( £ ) 


where I a 


AneKNLE 

"fmc 2 


( 22 ) 


Specifically, 

oo 

£ = phase acceleration of the electron, 

O 

£ = v = phase velocity of the electron, 
£ = phase of the electron, and 
| a | = optical field amplitude. 


• o 

We also know that £ = (k + k 0 )z - co implying that t, = L[(k + k 0 - k ]. Now 
we have a relationship between the phase £, and the phase velocity, v. 

The phase and phase velocity can be plotted on a phase space plot that shows 
the evolution of the electron trajectories in phase space. The path separating open 
and closed orbits in phase space is called the "separatrix", 

v s 2 = 2 | a | [1 + sin( ^ + 4> )] , (23) 

where 4> is the optical phase. The separatrix passes through the fixed points 
(-tc/ 2,0) and (3rc/2,0). The peak-to-peak height of the separatrix is 4 | a p /:2 , 
while the horizontal position is determined by the optical phase (Figures 8 and 9). 
The electrons are injected into the undulator and therefore start their phase space 
orbits with initial phase, £o> and initial phase velocity, v 0 . How well the electron 
distribution is arranged on the phase space picture is a direct indication of the 
beam quality. The poorer the beam quality, the broader the distribution will be. 
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Figure 9 . 

Schematic of the phase relationship between the electron 
trajectory and the light wave. 
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Figure 8 is a plot of phase space that shows graphically the relationship we have 
just developed. Figure 9 is a plot of the relative positions of the light, and the 
electron path as the phase space plot cycle begins. The numbers show possible 
relative positions between the light and the electron as they enter the undulator. 
The correspondence between the numbers in Figure 8 and Figure 9 are: 

Point 1 (relative phase position shown) is an unstable critical point. 

Point 2 (if it were at the origin in Fig. 9) corresponds to a stable 
critical point. 

Point 3 (if it were at the origin in Fig. 9) describes the same condition 
as point 1. 

The gain of the FEL comes from a bunching of the electrons as they interact 
with the light that is best depicted in phase space. As seen in Figure 10, the 
electrons in the left hand side of the separatrix have gained in phase velocity 
(absorbed energy from the light) and the electrons in the right hand side have lost 
in phase velocity (lost energy tc the light). 
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*** FEL Phase Space Evolution *** 

a =2 v =0 
o o 

0.001 

0 

0.1 

0 

-n/2 £ 3n/2 Ox 1 

Figure 10 . 

Phase space plot for electrons injected in the resonance condition. 



If the electrons are allowed to continue the interaction for a longer time, the motion 
will continue. Some electrons lose energy and some gain, so that there is no net 
energy transfer and the gain is zero. If the electrons are given a non-zero initial 
phase velocity, there is a difference in the net energy exchange and non-zero gain is 
achieved (Figure 11). It can be shown that the maximum gain is achieved at 
v 0 = 2.6. 
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Figure 11 . 

Phase space plot for electrons injected with a positive initial phase velocity. 


The proceeding discussion shows that the electron positions and velocities are 
critical to the performance of the FEL. If the initial electron velocities are random 
over a significant range, bunching of the electrons degrades and the gain is small. 
If non-uniform magnetic fields cause deflections from the undulator axis, the 
electron phase velocities are changed, bunching can be destroyed and the gain is 
small. Thus, the FEL undulator, which controls the magnetic field magnitude, 
length of interaction and uniformity of fields, is critical to the ultimate performance 
of the device. It is easy to see that optimization of undulator design is crucial to 
all applications of the Free Electron Laser. 
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B. UNDULATOR DESIGN 


Many Free Electron Lasers have been built incorporating various attributes. 
As in most endeavours, there is no one design that has all desirable traits 
incorporated. Thus, the researcher must decide which specific parameters are 
most important and design the optimum device to achieve it. There are two 
general classes of undulators available today. They are the linearly polarized and 
the circularly polarized classes. The linearly polarized class primarily uses 
permanent magnet structures while the circularly polarized undulator primarily 
uses current carrying coils to establish the magnetic fields. 

1. Linearly polarized undulator 

The permanent magnet structures produce a linearly polarized magnetic 
field of the form [5], 

B L = B m [ 0 , sin (k 0 z )cosh (k 0 y ) , cos {k 0 z )sinh (k 0 y ) ) . (24) 

The perfect electron trajectories are sinusoidal in the y-z plane with no motion in 
the x direction. Away from the undulator axis, the average transverse field 
strength increases in the x direction only, 

__ L B m \ k 2 x 2 

+ —-+ •" J • (25) 

This undulator provides focusing in the x direction, but not in the y direction. In 
practice, the experimenter uses other forms of magnetic focusing, such as a 
quadrupole lens or machining a slight parabolic curve in the normally flat magnet 
pole faces. The curved surface produces a smaller gap between the magnetic poles 
off axis, thus increasing the field strength off axis. The parabolic shape results in 
equal focusing in each radial dimension to maintain a circular electron beam cross 
section. The field from parabolic pole design with equal x-y focusing is 
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(26) 


The perfect electron trajectories on the axis are again sinusoidal in the y-z plane, 
but away from the axis the transverse field increases in both the x and y direction 
equally, 


fff = 


B k 


<2 


1 + 


k 2 (x 2 +y 2 ) 


(27) 


Notice that the magnitude of the average field seen by the electrons on axis is 
- B m H 2. The average transverse acceleration of the electrons is smaller due 
to the magnetic field reduction, and thus, the linear undulator gain is reduced. As 
we will see, this is not a problem with the helical undulator, but the construction 
techniques involved make the linear undulator much less complicated to build. 


2. Circularly Polarized Undulator 

The helical undulator with circular polarization has a field involving first- 
order Bessel functions, but near-axis fields can be approximated with error <1% , 
for k 0 r < 0.8 by [5,6] 


B X H = -B r 


1 + jk 0 2 (3x 2 + y 2 ) 
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, and 


B” = -B m k n 


1 + jk 2 (x 2 +y 2 ) 


[x cos( k 0 z ) + y sin(& 0 z )] 


(28) 


Electrons injected into these fields will travel in helical paths centered along the z 
axis. Comparing the average transverse magnetic fields, we find that 
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(29) 



1 + 


k 2 r 2 



and it is easy to see that off-axis fields are stronger closer to windings. This radial 
increase of the fields has a "focusing" effect on the electrons which results in a slow 
oscillation around the longitudinal axis. This motion is known as Betatron 
oscillations. The magnitude of the average field seen by the electrons on axis is 
= B m which results in higher transverse acceleration of the electrons. This in 
turn causes greater electron bunching, and thus, higher gain. 

Greater coupling and simplicity led us to choose the helical bifilar design. 
With this design in mind, another attribute we desire is a shorter wavelength light 
output. One of the ways to accomplish this is the electromagnetic compact FEL. 


C. THE COMPACT FEL 

The compact FEL design is one method of accomplishing shorter wavelength 
light output. Since we know that 

X = —( 1 + K 2 ) , (30) 

by shortening the length of our undulator period, we can shorten the light 
wavelength. Shorter wavelength light has many desirable properties such as 
atmospheric window propagation, and laboratory material interaction to name just 
a few. Shorter wavelength light also has higher energy so energy deposition on 
targets is greater. Thus, the compact bifilar helix is the design chosen to illustrate 
our modeling technique. 


1. Specifications 

Originally proposed by Roger Warren (LANL) [7], and considered as a 
tentative experiment for LANL and the Naval Postgraduate School, the following 
parameters were used: 
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Number of periods N = 10 , 

Undulator period X 0 = 0.9cm , 

Gap 2g = 4mm , and 
Lorentz factor y = 30.35 . 

Although this is the basis of comparison throughout this paper, all lengths were 
non-dimensionalized by dividing by X 0 for ease of transposition to other undulator 
designs. 

2. Winding schemes 

Using the dimensions for the proposed compact FEL, it is easy to see that 
the windings of the bifilar helix cannot be of a complicated nature. The small 
dimensions make the winding terminations at both ends of the undulator, and the 
method of power lead connection extremely important to the resultant magnetic 
field structure at the entrance to the undulator. We will see that the lead wires at 
the end of the undulator carry sufncient current to cause unwanted perturbations 
in the incoming electron beam. 

Three terminations are investigated in this paper: the wire, loop and 
staggered termination schemes [8]. In the case of wire termination, the bifilar 
helical wires are continued radially outward when they reach the end of the 
undulator as shown in Figure 12a. In the case of loop termination, a circle of wire 
is attached to the end of the cylinder formed by the helix with the windings 
attached at opposite positions on the circle. Thus, the current from one helical 
winding enters at one position on the circle, flows equally around each side of the 
circle, and then combines to flow in the opposite direction in the opposing helical 
winding as shown in Figure 12b. Staggered terminations are achieved by using 
wire or loop terminations as a current taper for the undulator removing a specified 
percentage of the current per termination depending on its position in front of an 
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extended undulator. This is visualized by extending the bifilar helix windings past 
the entrance to the undulator, and connecting them periodically by a termination 
which passes a percentage of the current in the helical windings. Thus, the current 
flowing in the helical windings is reduced stepwise by the number of terminations 
added prior to the entrance of the undulator. This results in a tapering of the 
magnetic field. Tapering of the magnetic field at the entrance of the undulator can 
also be achieved by flaring the windings which is accomplished by using a 
progressively larger radius for the extended bifilar windings as distance before the 
original undulator entrance increases. The Loop termination scheme with current 
tapering (staggering) vice flaring the windings results in the smallest magnetic 
field spike at the entrance of the undulator [8]. This paper investigates fields at 
the entrance of the undulator in the compact FEL. The effect of undulator exit 
termination is explored as well, and the fields found are used to determine the 
incoming electron trajectories. 
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III. MODELING OF UNDULATOR CHARACTERISTICS 

A. REPRESENTATION OF FIELDS 


1. Theory 

As discussed in the background section, this paper calculates ti e 
magnetic field at any point by using the Biot-Savart law. The magnetic field at 
any position can be expressed by, 


A# = MaTxT 
4k r 3 


where r = [ r 2 + r 2 + r 2 2 ] V2 . 

A& is the magnetic field contribution at the desired position by th current element. 
Po is the magnetic permeability constant. I is the current flowing in the current 
element, AL. The distance and direction of the desired position from the current 
element is expressed by T*. By integrating over the winding configuration of choice, 
the magnetic field at any particular position can be found. To find the current 
elements, All, we need to investigate the particular geometry of our winding 
configuration. 


a. Integration of Magnetic Field 

To develop the integration, we describe a very small piece of a 
current carrying wire in a particular geometry. As discussed previously, the 
methods of termination discussed will be; (i) Wire, (ii) Loop, and (iii) application of 


26 






wire and loop to the exit configuration. First, we need to discuss the geometry of 
the helical windings to set up the periodic magnetic field needed for undulator 
operation (Figure 13). 



r 

Figure 13 . 

Schematic representation of the bifilar helical path. 
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where g = radius of the undulator. 


Now use the definitions 

A L = (ALj. 2 + ALy + &L 2 ) V2 , and 
r z = - Xj ) 2 , r y 2 = (y -yj) 2 , r 2 = (z - Zjf , 

where ( x , y , z ) is the position of the electron, 

to substitute into (31) to sum up all the contributions from the current elements in 
the helical windings. If we did not have to worry about end effects, this method 
would give the desired solution to the interior magnetic fields of the undulator. 
Unfortunately, this is not an infinite length undulator and we next start 
discussing the termination methods. The simplest method of termination and 
simulating the leads connected to a power supply is the wire termination method 
(Figure 14). 




>1 = SL s2 = 

n z n z 

Figure 14 . 

Schematic representation of the wire termination method. 
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(2) Wire Termination. Referring to Figure 14, we can describe 
the position and angle of the current element mathematically by, 




a + e 


2ng 

n z 


_ ISL 

n 2 


y e = o , z e = o 


where e = 1 , 2 , • - • , 7 n z . 


The current elements have the form: 

A L x = , AL V = 0 , AL 2 = 0 . 

The contribution of the wire termination is summed by substituting into (31) as 
explained for the helical windings. The top and bottom wires need to be 
differentiated in their individual contribution. The current elements in both 
produce the same respective orientation but the distance from the current element 
to the electron is different according to the following: 

Topwire: r e = [ (x - Xg) 2 + y 2 + z 2 ] m , 

Bottomwire: r e = [ (x + x^) 2 + y 2 + z 2 ] V2 . 

As will be seen, this method of termination is by no means ideal. The other 
method of termination explored in this paper is loop termination (Figure 15). 
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Figure 15 . 

Schematic representation of the loop termination method. 


(3) Loop Termination. Referring to Figure 15, we can describe 
the position and angle of the current element mathematically by, 

x t = gcos0 e , y e = g sin0 e , z e = 0 , 

where 0 e = and e = 0, 1 , 2, ..., n z . 


The current elements have the form: 


„ 2jigsin0, 

0 e = 0 -» 180°: A L z = — -- , AL y = 

n z 


-2Kg COS0 e 


, AL z = 0 


„ -2^sin0. 

0, = 180 -> 360°: A L x = ---- , A L y 

n z 


2Kg COS0 e 


, A L z — 0 


The current elements are summed as before. The last termination investigated is 
applying the wire and loop terminations to the exit of the undulator. 
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(4) Exit Termination. To apply wire termination to the end of 
the undulator, the only difference in the previous argument is to reverse the sign of 
the current flowing in the current elements. To apply loop termination to the exit, 
and still connect the power leads, requires close examination of the geometry 
(Figure 16). 




X 



Figure 16 . 

Schematic representation of loop termination at the exit (up). 


As can be seen in Figure 16, the current in the loop exit termination is in the 
opposite direction. Also, replace all uses of z in the equations with z - N to get 
the correct distance to the electron. The wire termination used to attach the power 
lead to the loop is a modification of the previous geometry using z - N and the 
distance now becomes, 

r e = L (x -x e ) 2 + (z - N) 2 ] m . 

The wire termination used at the end of the second wire from the bifilar helix is 
the same as above except the current is reversed and the z distance becomes 
z - 0.999 N for physical separation of the leads. 
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The final variation studied is the exit termination with the 


wire connections going down (Figure 17). 



X 



Figure 17 . 

Schematic representation of loop termination at the exit (down). 


The only variation in the equations from the previous discussion is in the r x term. 
Replacing the minus sign with a plus sign achieves the necessary modification. 

(5) Coil Winding Tolerances. Though not pursued in this 
paper, there is another possible application of the integration method discussed 
above. The flexibility of the current element integration technique allows easy 
change of the winding configuration. This flexibility would apply just as easily to 
displacing the windings slightly from the ideal positions explored. This would 
result in the ability to compare simulations to determine coil winding tolerances. 
The usefullness of this would help in the construction of actual undulators as 
outlined in [9]. 
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b. Use of the K vector 


The magnetic field at any point in the undulator is found from (31). 
A more useful quantity in the derivation used for the electron trajectories is the 
value K. Since the definition for K is K = eB m X 0 /(2rcme), we see that the value 
of K is proportional to B. Thus, the representations of the magnetic field s^own 
will be using K to lead into the electron trajectory discussion. 

2. Results 

Now that we have developed the theory, we need to have some method of 
displaying the fields so we can get a better understanding of the stray field 
problem. As mentioned earlier, the Biot-Savart law was used to calculate the field 
contributions due to individual current elements. The field near the axis can be 
found analytically from (28) and agreed with the numerical calculations. Thus, the 
accuracy of the model is good, and can be applied to positions outside the 
undulator. The integration method also allows us flexibility in the winding 
configuration used for the integration. Thus, we can model the stray fields at the 
entrance to the undulator with the various termination schemes already discussed. 
The integration technique is first used to find the magnitude of the K vector along 
the axis for comparison between winding configurations, and as a comparison to 
the original work by Fajans [8]. In all simulations, the undulator entrance is fixed 
at z = 0. If tapering is used, it is added prior to this point making the undulator 
longer. 


a. K Components 

The first simulation is the K values for wire entrance and exit 
terminations. This will be referred to as wire/wire termination. This is the 
simplest of the termination techniques explored. The method of display is a dual 
plot of the x and y components of the K vector along the z axis. The z scale shows 
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the tapering used by a circle whose radius is proportional to the current centered 
on the undulator periods (Figure 18). 

There is a slight asymmetry to the magnitudes of the two 
components as seen in the maximum and minimum values displayed. Notice that 
the y component of the K vector has a significant value for a considerable dis 1 ace 
before z = 0. This is expected as the wires oriented for the termination contribute 
a field on axis in the y direction only. There is no large field gradient or spike at 
the entrance as predicted by Fajans [8]. The distance over which the y component 
can affect the entering electron trajectory may result in the same deflection as a 
large field spike. 


** K magnitude for wire/wire terminations ** 

„ Kxmax=l.04 Kvmax=l.28 

N=10 Kxmin=-1.05 Kymin=—1.00 

Hi 

mm 
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Figure 18 . 

Combined plot of the K components along the z axis with no taper. 
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Also, there is an abrupt transition in the magnitude of the fields, 
especially in the x plane, right at the entrance to the undulator. To reduce the 
abrupt change, we try tapering the undulator. The first example of tapering uses a 
1 period taper shown in Figure 19. 


** K magnitude for wire/wire terminations ** 

, - Kxmax=l.04 Kymax=l.21 

:1 ° Kxmin=-1.05 Kymin=-1.00 



-N/2 Z N/2 

Figure 19 . 

Combined plot of the K components along the z axis with 1 period taper. 


As is seen, tapering ameliorates the abrupt field changes at the entrance. 
Unfortunately, the long lead of the y component is still evident. Also, even though 
tapering has been applied, there is a significant asymmetry still evident just inside 
the entrance. 
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This effect does not go away, even if the tapering is increased as 
shown in Figure 20. Like before, the long tail and asymmetric spike just inside the 
entrance is seen. The asymmetry of this termination scheme leads us to the 
suspicion that this winding design might not prove to be optimum. 



Figure 20 . 

Combined plot of the K components along the z axis with 2 period taper. 


The next step is to try loop termination to see if better results are 
achieved. If we apply a loop entrance and a wire exit (loop/wire) termination, we 
arrive at the values displayed in Figure 21. We see immediately that the range 
over which the y component is a single contributor is considerably shorter. The y 
component also starts with a positive value vice a negative one. This shortening 
and reversal is accompanied with a better symmetry seen throughout the displayed 
range as the maximum and minimum values show. 
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** K magnitude for loop/wire terminations ** 
Kxmax«=l. 02 Kymax=l. 03 

N=10 Kxmin=-1.03 Kymin=-1.04 
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Figure 21 . 

Combined plot of the K components along the z axis with no taper. 


The improved symmetry should get better as we add tapering and it 
does. Figure 22 shows the combined components with a 2 period taper. The short 
range that the y component acts unopposed outside the entrance, and improved 
symmetry should prove to be a better design for the injection of electrons which 
will be explored later. 
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K magnitude for loop/wire terminations 
Kxmax=l.02 Kvmax=l.03 

Kxmin=-1.03 Kymin=-1.02 
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Figure 22 . 

Combined plot of the K components along the z axis with 2 period taper. 


The improved symmetry seen for the loop termination can be taken 
one step further by using it on the exit termination. This loop/loop termination 
would simulate an undulator with loop termination at both ends, and the power 
leads connected to the power supply via coaxial cable that give no magnetic field. 
We expect that this modification should have slightly better symmetry, since the 
dimensions of the compact FEL bring the exit termination close enough to have a 
noticeable effect on the entrance fields. If loop exit termination is added to the 
previous simulation, mixed results are achieved as shown in Figure 23. The range 
between minimum and maximum values for the components is reduced and the 
minimums and maximums are more constant than the previous case. 
Unfortunately, there is no clear cut indication that this termination technique 
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produces noticeable changes. Thus, we can only guess that this method would 
provide a better entrance field for incoming electrons. Since we have verified 
correlation to previous work [8], and established a guess as to which one should 
offer the most symmetrical magnetic fields, we need to find a better way of 
representing the fields. This section investigated the field amplitudes along the 
axis. The electrons will be following a helical trajectory through the undulator, and 
therefore, will experience fields off-axis in their travel through the undulator. 
Thus, modeling of off-axis fields needs to be explored. 


** K magnitude for loop/loop terminations(up) ** 

„ __ Kxmax=l.00 Kymax=0.99 

N-10 Kxmin=-1.01 Kymin=-1.01 



-N/2 * N/2 

Figure 23 . 

Combined plot of the K components along the z axis with 2 period taper. 
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b. Closeup of Entrance K Field 

As mentioned previously, the fields off-axis are important to the 
calculation of the electron trajectories. A model can be envisioned that would 
contain the space vectors representing the magnetic field at any point in the 
undulator. This is the definition of a vector field, and while the subsequent 
electron beam trajectory calculations would be relatively simple, it is easy to see 
that any attempt to give accuracy to the vector field representation results in a 
computational nightmare. The sheer number of points needed for a three 
dimensional representation would overwhelm most computers. Thus, the model 
needs to be reduced in scope. The most obvious choice is to limit the boundaries of 
calculation inside the undulator which will reduce the number of points in the 
vector field that need to be calculated. Since the dimensions of the electron beam 
and the limits of off-axis motion are small, this is a reasonable simplification. The 
boundary could be as small as a fourth of the radius of the undulator and the 
number of points would then be reasonable for present-day computers to calculate 
and display. This method would be the best to pursue if a real beam of electrons is 
trying to be simulated. As mentioned in the background section, the density of the 
beam is on the order of 10 12 electrons /cm 3 . It is not practical to calculate the 
magnetic field for individual electrons when this many are involved. However, the 
model in this paper does not deal with a great number of electrons, only a small 
representative fraction. The calculations can be greatly speeded up by only 
calculating the magnetic field at the electron positions. This drastically reduces 
the number of calculations needed for accurate representation. The magnetic fields 
at the electron positions are used to calculate the force. This enables us to find the 
accelerations, and subsequent integrating finds the new positions. 

Now that we have established that calculating a complete vector 
field is too hard, it would be even worse trying to represent it in some intelligible 
way. Therefore, we will pursue a planar representation of the vector field. The 
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integration technique was applied to the end of the undulator to see the stray field 
effect. The structure of the fields were illustrated by shade gradations 
corresponding to the value of the K vectors. The magnetic field values get very 
high close to the windings as one would expect. Therefore, to see the fine structure 
of the fields, positions only out to half the radius of the undulator are shown. The 
magnitude plot on the left has the windings on the outside of the undulator 
overlaid as dark black lines to help provide viewpoint recognition. To emphasize 
the need for tapering, simple wire termination is shown with positions from one 
period outside the undulator to one period inside the undulator as shown in 
Figures 24a and 24b. The white lines are contour lines of constant magnitude. 
Notice the skewed nature of the fields from the contour lines. As was conjectured 
previously, there are large asymmetrical fields at the entrance, and undesirable 
perturbations would distort the incoming electron trajectories. 
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*** K values vs. position 
No taper X o =0.009m g=0.2222 



IK(0,0,0)| = 0.4691 


Figure 24a . 

Graphical representation of the K field composition in the x plane 
from 1 period outside to 1 period inside the undulator. 

The design is wire/wire termination. 
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*** K values vs. position *** 
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Figure 24b . 

Graphical representation of the K field composition in the y plane 
from 1 period outside to 1 period inside the undulator. 

The design is wire/wire termination. 
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If we rerun the same scenario with a two period taper, we see by the 
contour lines and the maximum and minimum values that some of the 
asymmetrical nature has been reduced as shown in Figures 25a and 25b. Even 
though tapering has improved the picture, the overall symmetry of the wire 
termination scheme still leaves a lot to be desired. The next step is to apply loop 
termination to our model and see if our previous conclusion of better symmetry 
applies to off-axis fields. 

Overall symmetry is seen when we apply loop/wire termination to 
our simulation as shown in Figures 26a and 26b. Even with no taper applied, the 
contour lines and the maximum and minimum values show a much more organized 
nature than before. This trend should get better if we apply tapering to our model, 
and it does in Figures 27a and 27b. The final simulation will be for loop/loop 
termination with no taper for comparison to the loop/wire termination (Figures 28a 
and 28b). 
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Figure 25a . 

Graphical representation of the K field composition in the x plane 
from 1 period outside to 1 period inside the undulator. 

The design is wire/wire termination. 


|K(0,0,0)| = 0.7414 
Two period taper 


g/2 —g/2 x g/2-g/2 x g/2 


45 













0.40 


| K(0,0,0) | = 0.7414 
Two period taper 


Figure 25b . 

Graphical representation of the K field composition in the y plane 
from 1 period outside to 1 period inside the undulator. 

The design is wire/wire termination. 
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*** k values vs. pc 
:aper X =0.009m 


Figure 26a . 

Graphical representation of the K field composition in the x plane 
from 1 period outside to 1 period inside the undulator. 

The design is loop/wire termination. 
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*** K values vs. position *** 
:aper X =0.0 09m g=0.2 


Figure 26b . 

Graphical representation of the K field composition in the y plane 
from 1 period outside to 1 period inside the undulator. 

The design is loop/wire termination. 
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|K(0,0,0)| = 0.8508 
Two period taper 


Figure 27a . 

Graphical representation of the K field composition in the x plane 
from 1 period outside to 1 period inside the undulator. 

The design is loop/wire termination. 












*** K values vs. position *** 
Tapered X =0.009m g=0.2222 



| K(0,0,0)| = 0.8508 
Two period taper 


Figure 27b . 

Graphical representation of the K field composition in the y plane 
from 1 period outside to 1 period inside the undulator. 

The design is loop/wire termination. 
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*** K values vs. position *** 

No taper X Q =0.009m g=0.2222 



-g/2 x g/2-g/2 x g/2 

1K(0,0,0) | = 0.5458 


Figure 28a . 

Graphical representation of the K field composition in the x plane 
from 1 period outside to 1 period inside the undulator. 

The design is loop/loop termination. 
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Figure 28b . 

Graphical representation of the K field composition in the y plane 
from 1 period outside to 1 period inside the undulator. 

The design is loop/loop termination. 
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Unfortunately, the exit termination representation does not give a 
clear improvement of the stray field structure. The resolution of optimum design 
will have to wait for a different representation of the fields. To pursue this avenue, 
there is a possibility that the view shown is too narrow. Thus, we shall expand the 
displayed K field to show an enlarged view. 

c. K Field Over Five Periods Inside and Out 

We next apply the terminations to a display that shows positions 
along the z axis from 5 periods outside the undulator to the middle of our 
undulator. Again, shade gradations show the detailed structure of the stray fields 
at work on the entrance to the undulator as shown in Figure 29. The structure 
seen is quite complicated and several trends are evident. The first is the contour 
lines along the axis of the undulator. The contour lines show that the fields are 
not uniform along the axis where they should be the most symmetrical. Also, the 
light area right at the entrance to the undulator indicates the strongest magnetic 
field is present right at the entrance. This abrupt field change may prove 
detrimental to the electron trajectories investigated shortly. If the contour lines 
close to the coil positions are studied, the x plane reveals an interesting feature. 
The contour lines have a magnitude sinusoidal with the z direction. The 
magnitude of the contours gets increasingly symmetrical as the undulator is 
traversed indicating the termination wires affect the magnetic field structure all 
the way to the center of the undulator. This effect is not evident in the y plane. 
We could conclude that electron trajectories would be significantly affected by this 
type of design, especially in the x plane. Since the stray field structure looked 
more symmetrical for loop termination, we should see if this applies to our 
expanded view of the undulator fields. 
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** Dndulator Magnetic Fields ** 
N=10 gA o =0.2222 K<0,0,N/2)=1.00 
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Figure 29 . 

Graphical representation of the K field composition in the x and y plane 
for wire/wire termination. The displayed positions are from 5 periods 
outside the undulator to 5 periods inside the undulator. 


If loop/wire termination is applied to our model undulator, we can 
see major changes in the representation of the fields as shown in Figure 30. The 
overall symmetry of the structure is better, especially in the x plane. The 
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sinusoidal variance of the fields close to the windings now appears to be reversed, 
but at a significantly lower amplitude. The abrupt change in field intensity at the 
entrance is also reduced in intensity. The close spacing of the contour lines at the 
entrance signify a significant field gradient right at the entrance that seems to 
extend for a much shorter distance than wire/wire termination. An interesting 
feature is the significant asymmetry at about 1 undulator period inside the 
undulator. This is evident in both the x and y plane. The conclusion taken from 
this representation is that this design provides better magnetic field symmetry and 
should have a less detrimental effect on the incoming electron trajectories. The 
final design considered is for loop/loop termination as shown in Figure 31. 
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** K values entering undulator ** 
N=10 g/X o =0.2222 K(0,0,N/2)=1.00 
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Figure 30 . 

Graphical representation of the K field composition in the x and y plane 
for loop/wire termination. The displayed positions are from 5 periods 
outside the undulator to 5 periods inside the undulator. 


Looking at the plot, we still see the ambiguity we had before with 
this design. Since the construction of the undulator is close to loop/wire 
configuration, it is not surprising that we see the basic symmetry throughout the 
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displayed fields. The major difference seen is a sharp reduction in the sinusoidal 
variations of the contour lines close to the win din gs. The increased symmetry seen 
should prove to be the best design of the termination schemes explored. 
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** K values entering undulator ** 
N=10 g/X o «0.2222 K<0,0,N/2)=1.00 



—N/2 Z N/2 


Figure 31 . 

Graphical representation of the K field composition in the x and y plane 
for loop/loop termination. The displayed positions are from 5 periods 
outside the undulator to 5 periods inside the undulator. 
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B. REPRESENTATION OF TRAJECTORIES 


1. Theory 

As mentioned previously, the magnetic field is found a; the electron 
position only when calculating the electron motion. Using the magnetic field, the 
force is found using the Lorentz Force equation (5). We neglect the presence of 
light ( E r = = 0 ), and use the calculated magnetic field to integrate the motion. 
We know that 


-j-(yj?) = ——[ (? x £? ] and if y = 0 
at m 


( 6 ) 


=> ft. = -£ z B y ,%B X -p x B z ,0] , 

ym 

P 2 = ~[ £ z B y - £ y B x ] . (32) 

Recall from (17) that 

P x - — <c 1 , (33) 

Y 

so that p z = 0. Then, using P z = 1, we have the approximate z motion, 
z = z 0 +tL where x = ct/L. Using p iy c p z , we also find that 


which leads us to 




,-B x ] 


(34) 


x_ 

c 



c 



(35) 


Using dt - Ld ilc , we find that 


oo 

X 


L 2 eB y 

ymc 


00 

and y 


L 2 eB x 

ymc 


(36) 


o 

where ( ) = d{ )/dx. Finally, voth x/X Q -* x ,y/X 0 —> y , and z/X 0 -->z, we 
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simplify these expressions to 


00 2j zN 2 K„ 

x = -- 

Y 


where 


OO 

y 


2 kN 2 K x 

y 


K. 




2 Time 


z = z 0 +Nx , (37) 


The Euler-Cromer integration method will be used to update the electron position 
as it travels through the fields represented in our simulated undulator [10]. 


2. Example of Ideal Injection 

The previous discussion explores the technique used to determine the 
electron trajectories. This .'action describes ideal injection where the electron is 
moving tangent to the motion it would describe if it were inside of an infinite 
undulator. In our case, this motion is a helix centered on the z axis. Good 
injection is desirable in experiments because it establishes a smooth transition into 
the undulator and increases the gain of the FEL. 

The exact trajectories in a helical undulator are 



-cos(o ) 0 t) ,-sin( 0 ) o O 

Y Y 



(38) 


where 


0 ) o = P 0 k 0 c = and P 0 2 = 1 - 1 + K 


h 


Y 2 


By integrating, we cam get the exact equation for the electron position at any time 
t, 


r K\ n 


Kk„ 


ot = —^-sin(o) 0 f) , 

L 2x7p 0 2 xtP 0 


cos(co 0 t) , P c ct 


(39) 


At time t = 0, the trajectory describes the maximum angle with the z axis in the x- 
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z plane. Therefore, to describe the angle that is tangent to the helix, we define the 
angle 8, where it can be shown that, 

e x = — p,(0) = , 6^=0 if (3,(0) = p o = 1 , 

u>z Y 

and 


Ax = 0 


KX 0 

2tcyPo 


(40) 


By substituting the parameters of our model undulator with K = 1, y = 30.35 
and X a = 0.9cm. into (40), we find that the perfect injection angle is G z = 0.033 
radians, and Ay = 0.0053 (non-dimensional units). 


3. Results 

Now that we have developed the theory of electron trajectory calculation, 
we need to develop some way of displaying the resulting information. To add some 
continuity, the electron trajectories will be calculated and overlaid on the same K 
field shown previously in Figures 41, 43, and 45. Although simple in principle, 
there are a couple of complicating features that need to be addressed. The first is 
the electrons’ initial position and angle inside the beam. The second is finding 
some way of estimating the injection parameters so that optimization of undulator 
designs can be discerned. The easiest way to accomplish this is to determine a 
"characteristic" angle and position offset that the stray magnetic fields will bend 
the individual electrons through. 

The individual electrons inside the simulated electron beam are given an 
initial position and angle consistent with a characteristic beam quality factor called 
"emittance". The definition for emittance is e* = 2 k x 8, and Zy = 2 k y Q y , and 
it is assumed in this paper that e, = Zy. The quantity x is the rms initial 
position spread of electrons along the x direction and 0, is the rms initial angular 
spread of elect* *ns from the axis of the undulator along the x direction. Either the 


61 






rms position spread, x, or the rms angular spread, 0 X , can be changed by 
external focusing fields prior to entrance into the undulator, but their product, £*, 
is fixed. Thus, arbitrary values for X and 6 Z are chosen such that X = 2 r e , 
where r e is the electron beam radius, and 0* =4 r e /N. The value used in this 
paper is r e = 0.1 g where g = 2mm. 

The electrons are initially started at 5 undulator periods away from the 
beginning of the first taper. The electrons are injected at a position and angle such 
that the resulting trajectories execute helical motion centered around the z axis. 
Since there are a large number of designs, we need to find some method of 
determining the optimum tapering scheme for each termination. The initial angles 
and positions are estimated by a calculated characteristic deflection angle and 
position offset for each design. By finding the minimum characteristic deflection 
angles and position offsets, the optimum taper for each end termination is 
determined. Using (37), we have 

°° 2nN 2 K v dv x , 06 2kN 2 K x dv Y 

x = -— = —— and y --- —~ 

y dx y ax 

With dx = dz/N, it can be shown that on the undulator axis, 

- ^C»^ <o ' o ' z)cb • 

where t is the number of tapered periods. The initial transverse velocity gives the 
initial transverse angle, so that 

dv x = = e N . (41) 

dx dz 

Therefore, the characteristic angle of an electron traveling into the undulator is 

e * = Y^^ Ky{0,0,z)dz and Qy = -f£^ K ' m ’ z)dz ' (42) 

To find the characteristic position offset, we continue the previous derivation and 
apply (41) to (42), and integrate dx over the length of the trajectory. The resulting 
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offset is given by 


*= = Yl<w)*'C»^ <0 ' 0 ' z)dz 31,(1 

^ - - 7 ■Lu dz 'C,^ wfi - z)dz • <43) 

By translating (42) and (43) into a computer program, the different 
terminations and tapers can be evaluated in tabular form shown in Tables 1 to 3. 
The type of tapering explored was either linear or a smooth curve fit to the 
successive current reductions. If more than a three period taper is used, there is a 
possibility of designing the taper such that the stepwise current reductions do not 
follow a straight ramping function. The current reductions can be tailored to follow 
a curved function. The curve fit was tried with both a tangent function and a 
smooth curve (determined by hand) representing the current reductions. 

It is important to note several points about the characteristic angles and 
positions determined. The first is that the K values are confined along the z axis. 
Since the actual electron trajectories do not exactly follow the z axis, the deflection 
angles and position offsets found will not exactly match the offsets required for 
injection. However, the trends found in the tables are a good indicator of optimum 
tapering schemes for each termination design. The second point is the integration 
length of the y component. Referring back to Figures 18 to 23, notice the difference 
between the component K values along the z axis. The K x component ends on an 
even period at z = 5. This is not true for the K y component. The integration 
distance for the y component is taken to be one-quarter period less to account for 
the alternating values seen by the incoming electron. Therefore, the modified 
characteristic angle becomes 

e * * Yl<w)^ (0 - 0 > 2)d2 “rf e > = - yO ^ (0 ’°’ z)d2 • (44) 
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The new characteristic position offset becomes 

9 _ 4.75 z' 

A -y = -\li M) dz ’l.M) K * m ’ z)dz > (45) 

where K = 1.0 and y = 30.35. The units used for the tables will be radians for 
angles and non-dimensional lengths for offsets. 

The most notable feature in Table 1 is that more tapering does not 
achieve the best entrance conditions [8,9]. We can see that no tapering achieves 
the minimum characteristic values. Also, the characteristic angles are fairly 
insensitive to tapering. This is especially true for the y values. Recall that 
tapering did little to assuage the asymmetry at the undulator entrance. Although 
not easy to discern, it appears that if tapering is applied, the best method is a 
linear taper. This is another surprise as a smooth current taper should produce 
the smoothest magnetic field transition. Therefore, the optimum winding 
configuration for wire/wire termination is with no tapering. At the bottom of Table 
1 are the absolute values of the ideal injection angle and position offset at z = -5 
for the perfect trajectories shown in (38) and (39). The absolute values are used 
because the simulation and ideal polarizations are not matched. 
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Characteristic values for various tapering schemes with wire/wire termination. 


Table 1 . 

Tapering scheme 

o. 

Q y 

Ax 

Ay 

No Taper 

0.062 

-0.003 

-0.109 

-0.188 

1 period taper 
(linear) 

0.063 

-0.003 

-0.122 

-0.207 

2 period taper 
(linear) 

0.063 

-0.003 

-0.135 

-0.225 

3 period taper 
(linear) 

0.063 

-0.003 

-0.146 

-0.244 

3 period taper 
(tangent function) 

0.063 

-0.003 

-0.153 

-0.243 

3 period taper 
(smooth curve) 

0.063 

-0.003 

-0.150 

-0.243 

4 period taper 
(linear) 

0.064 

-0.003 

-0.156 

-0.262 

4 period taper 
(tangent function) 

0.064 

-0.003 

-0.171 

-0.261 

4 period taper 
(smooth curve) 

0.064 

-0.003 

-0.166 

-0.261 

5 period taper 
(linear) 

0.064 

-0.003 

-0.165 

-0.280 

5 period taper 
(tangent function) 

0.064 

-0.003 

-0.1° 1 

-0.279 

5 period taper 
(smooth curve) 

0.064 

-0.003 

-0.178 

-0.280 

1 Ideal values 1 

0.033 

0 

0 

0.170 


































































Next try the characteristic test on loop/wire termination to find if our 
original assumptions are true as shown in Table 2. 
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Characteristic values for various tapering schemes with loop/wire termination. 


_ Table 2 _._ 

Tapering scheme 0* 0 y Ax Ay 

No Taper 0.056 -0.003 0.108 -0.184 

1 period taper 

(linear) 0.057 -0.003 0.117 -0.203 

2 period taper 

(linear) 0.059 -0.003 0.124 -0.221 

3 period taper 

(linear) 0.059 -0.003 0.132 -0.239 

3 period taper 

(tangent function) 0.059 -0.003 0.132 -0.239 

3 period taper 

(smooth curve) 0.059 -0.003 0.132 -0.239 

4 period taper 

(linear) 0.060 -0.003 0.139 -0.258 

4 period taper 

(tangent function) 0.060 -0.003 0.140 -0.257 

4 period taper 

(smooth curve) 0.060 -0.003 0.140 -0.257 

5 period taper 

(linear) 0.061 -0.003 0.146 -0.276 

5 period taper 

(tangent function) 0.061 -0.003 0.147 -0.276 

5 period taper 

(smooth curve) 0.061 -0.003 0.147 -0.276 

i Ideal values I 0.033 0 0 0.170 









































































The most notable feature in Table 2 is that despite the enhanced symmetry in the 
fields noted with loop termination, no tapering is still indicated as optimum. 
Because of the increased sy mm etry, all values can be seen to be less than their 
wire/wire counterparts. The characteristic position offset in the x direction has 
reversed sign. As in wire/wire termination, if tapering is applied, the linear taper 
is slightly better. It is surprising to note that despite the increased symmetry, the 
values obtained are not significantly different from the wire/wire values. This 
implies that the electron trajectories will show an improvement in imposed input 
conditions, but not to any drastic extent. We will see if this is true when the 
simulation is carried out. Again, at the bottom of Table 2, the ideal injection angle 
and position offset are provided for reference. 

The last characteristics determined are for loop/loop termination as 
shown in Table 3. Table 3 shows a repeat of all the trends seen for loop/wire 
termination. As anticipated, the increased symmetry has resulted in smaller 
characteristic values for all methods explored. The inference here is that loop/loop 
termination will provide the smallest position and angle offsets required for 
optimum electron injection. The characteristic values found are close to the ideal 
values shown at the bottom of the Table. All terminations indicate that the 
optimum tapering scheme is no taper at all. We now need to apply this knowledge 
to an electron trajectory simulation. 
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Characteristic values for various tapering schemes with loop/loop termination. 


Table 3 . 

Tapering scheme 

0, 

0, 

Ax 

Ay 

No Taper 

0.037 

-0.003 

0.042 

-0.180 

1 period taper 
(linear) 

0.037 

-0.003 

0.041 

-0.198 

2 period taper 
(linear) 

0.037 

-0.003 

0.040 

-0.216 

3 period taper 
(linear) 

0.037 

-0.003 

0.038 

-0.234 

3 period taper 
(tangent function) 

0.037 

-0.003 

0.038 

-0.234 

3 period taper 
(smooth curve) 

0.037 

-0.003 

0.038 

-0.234 

4 period taper 
(linear) 

0.037 

-0.003 

0.036 

-0.252 

4 period taper 
(tangent function) 

0.037 

-0.003 

0.037 

-0.252 

4 period taper 
(smooth curve) 

0.037 

-0.003 

0.037 

-0.252 

5 period taper 
(linear) 

0.037 

-0.003 

0.034 

-0.270 

5 period taper 
(tangent function) 

0.037 

-0.003 

0.035 

-0.269 

5 period taper 
(smooth curve) 

0.037 

-0.003 

0.035 

-0.270 

1 Ideal values 1 

0.033 

0 

0 

0.170 
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Now that the optimum tapering scheme for each termination has been 
determined, we simulate the electron beam trajectories. The electrons are started 
at the characteristic deflection angle and position offset indicated for the optimum 
tapering configuration. Several simulations are run to minimize the input 
conditions necessary to achieve helical motion centered down the z axis. The first 
example is wire/wire termination shown in Figure 45. 

Figure 45 shows the superposition of the electron beam trajectories on the 
same design as in Figure 42 which has wire termination used at the entrance and 
exit of the undulator. The beam must enter in the correct position and angle in 
order to get the beam to stay in the undulator. After some optimization, the input 
position and angle parameters, starting back at z 0 = -5 away from the undulator 
entrance, are found to be x 0 = 0.45g , y a = -0.8 <g , Q x = 0 and Q y = 0.035 
where g = 2mm. When additional tapering is used, the input parameters become 
even more extreme. If no offsets are used, the asymmetry in the entrance fields 
sends the beam into the mde of the undulator at N - 4 periods. The stray fields 
resulting from wire termination at both ends are not suitable for practical FEL 
application. Loop termination provides a more symmetrical field composition, and 
should give a smaller deflection to the incoming electron beam. 
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If loop termination is used at the entrance of the undulator and wire 
termination is used at the exit, the result is Figure 46. The input angle 
parameters for this design, starting at z 0 = -5, become Q z = -0.027 and 
Q y = 0.035. The loop termination reduces the large magnitude field reversal in the 
x plane experienced by the electrons traversing the undulator entrance and thus, 
the resultant trajectories are much more like those shown in the y plane (electron 
beam deflected by a single angle at the entrance to undulator). Thus, only the 
angle offsets are given. The electron beam entrance requirements are reduced, 
especially in the x plane as predicted by Fajans [8]. Unfortunately, the y plane 
still suffers from about the same extreme position and angle offsets. It is possible 
that, because of the compact FEL dimensions, the exit leads are imposing an 
unexpectedly significant effect on the entrance stray fields. If loop termination 
could be applied to the exit leads, and the wires taken from the undulator in a 
coaxial cable, the effect of the wire exit termination would be reduced, and may 
ameliorate the entrance parameters. 
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** Electron trajectories using K ** 
N=10 ^30.35 g/X o =0.2222 #part.=25 


Figure 46 . 

Electron trajectories overlaid on K field background 
for loop/wire termination, no taper. 


Loop termination is applied to the undulator entrance and exit leads 
shown in Figure 47. T t is found that the minimum angle and position offsets are 
achieved when no tapering is used. The input parameter angles, starting from 
z o = —5, become 0 T = -0.012 and 0 V = 0.035. The result is that the input 
























electron beam parameters are nearly at the ideal injection angle of 0 ; 
calculated previously. This is obviously the best winding design. 


























IV. CONCLUSION OF OPTIMUM DESIGN 


The aim of this paper was to do a computer simulation of the fields and 
electron trajectories associated with an FEL. The compact FEL design used was 
originally proposed as a method of achieving shorter wavelength light output from 
an FEL. The simplicity of design lent itself to easy modeling and because of the 
small dimension:;, was assumed to accentuate the stray fields at the end of the 
undulator. These stray fields have not been extensively investigated as they are a 
small effect in more common FELs. This turned out to not be the case with the 
compact FEL. This determination was reached with many different types of 
representations of the magnetic fields associated with the device. The minimum 
angle and position offset imposed on the electron beam for smooth entry down the 
center of the undulator was used as the determining factor in the selection of the 
optimum design modification. 

All the design modifications considered show that stray fields at the ends of 
bifilar helical undulators with realistic wire leads cause serious deflections of the 
entering electron beam. Even with a termination that achieves improved 
symmetry in the magnetic field structure, loop termination, the discontinuity in the 
fields imposed by the termination itself with any kind of current tapering results in 
significant perturbations to the electron trajectories. The optimum design of those 
presented is the loop termination design at the undulator entrance and exit with 
no taper. This design achieves the smallest position and angle offset requirements 
on the input electron beam, and would be the least complicated to construct. 


75 






LIST OF REFERENCES 


[1] Pines, David, editor, et. al., Reviews of Modern Physics, Vol. 59(3), Part II, 
July 1987, p s33-s35. 

[2] Siegman, Anthony, Lasers, University Science books, 1986, p 2. 

[3] Colson, William, Free Electron Laser Theory , PhD thesis at Stanford 
University, 1977. 

[4] Hecht, Eugene, Optics, Addison-Wesley Publishing Co., 1987, p 272. 

[5] Colson, W. B., Pellegrini, C. and Renieri, A., editors, Free Electron Laser 
Handbook, Chapter 3, p 6,7. 

[6] Blewett, John P. and Chasman, R., Orbits and Fields in the Helical Wiggler, 
Journal of Applied Physics, Vol. 48(7), July 1977, p 2698. 

[7] Warren, Roger W. (LANL), private communication September 1989. 

[8] Fajans, J., Journal of Applied Physics, Vol. 55(1), January 1984, p 43-50. 

[9] Elias, L. R. and Madey, J. M., Superconducting helically wound magnet for the 
free-electron laser, Rev. Sri. Instrum., Vol. 50(11), November 1979, p 1335- 
1340. 

[10] Gould, H. and Tobehnik, J., An Introduction to Computer Simulation 
Methods Applications to Physical Systems Part 1, Addison-Wesley 
Publishing Co., 1988, p 38. 


76 







INITIAL DISTRIBUTION LIST 


1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, Virginia 22304-6145 

2. Library, Code 52 2 

Naval Postgraduate School 

Monterey, California 93943-5002 

3. Professor William B. Colson, Code PH/Cw 9 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93943-5000 

4. Professor Fred Ramon Buskirk, Code PH/Bs 1 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93943-5000 

5. Professor K. E. Woehler Code PH/Wh 1 

Chairman, Department of Physics 

Naval Postgraduate School 
Monterey, California 93943-5000 


77 





